Hierarchical Gompertzian growth maps with application in astrophysics 
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The Gompertz model describes the growth in time of the size of significant quantities associated to a large 
number of systems, taking into account nonlinearity features by a linear equation satisfied by a nonlinear func- 
tion of the size. Following this scheme, we introduce a class of hierarchical maps which describe discrete 
sequences of intermediate characteristic scales. We find the general solutions of the maps, which account for 
a rich set of possible phenomena. Eventually, we provide an important application, by showing that a map 
belonging to the class so introduced generates all the observed astrophysical length and mass scales. 
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PACS numbers: 89.75-k, 05.20.-y, 05.10.Gg 

The laws of growth of many systems, and the deep origin of 
their characteristic scales of, e. g., length, mass, energy, or nu- 
merosity, are intensively investigated in many branches of sci- 
ence, such as biomedicine, economy, and population dynam- 
ics. The Gompertz model was originally introduced in 1825 
by B. Gompertz lH as a model of human mortality: Gom- 
pertz found empirically the fitting distribution of human age 
for a given community. From the first half of the twentieth 
century the Gompertz model, and the associated Gompertz 
equation, have become a frequently used tool to account for 
mechanisms of growth 101 ■ The form of the Gompertz equa- 
tion is: 



at z 



(1) 



where z describes the "size" (not necessarily a spatial size) of 
some quantity characterizing the system, /3 and a denote two 
positive constants with the dimensions of the inverse of time, 
and 2 is a constant with the same dimensions of z. Eq. ([T]i can 
be recast as: 



rf(lns) 
dt 
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(2) 



where s{t) = z{t)/zoo, and Zoo = zexp{/3/a). The Gom- 
pertz equation is then characterized by four parameters (all 
dependent on the specific system): a, /3, z^o, and the initial 
condition ("scale") z(0) = zq. Its solution is: 



zit) = 2:00 cxp [(In 7) 



(3) 



where 7 = zq/zoo- It is immediately verified that this solu- 
tion always approaches monotonically in time Zqo- Depend- 
ing on the conditions zq < z^o = zcxp/3/a (7 < 1), or 
Zq > Zoo = zcxp/3/a (7 > 1), the system monotonically 
grows or monotonically decreases, respectively, from the di- 
mension Zq to the dimension z^o, approaching the asymptotic 
value Zoo with the characteristic time a^^. In TABLE I we list 
the main symbols which have been exploited, or which will be 
exploited in the following, with their meanings. 



Symbol 


Meaning 


a, 13 


Parameters of the Gompertz eq. 


z(t) 


Size at time t 


zo 


Initial size 


Zoo 


Asymptotic (or reference) size 


7 


Ratio ZQ / Zoo 


s{t) 


Relative size z{t)/zoo 




n-th relative size z„/zoo in the discrete map 


yn 


InSn 


a 


Parameter of the discrete map 


Rn 


n-th astrophysical length scale 


R 


Max. astroph. length scale (Radius of the obs. univ.) 


A 


Min. length scale i?o (radius of a nucleon) 


M„ 


n-th astrophysical Mass scale 


M 


Max. astroph. mass scale (Mass of the obs. univ.) 


m 


Min. mass scale A/o (mass of a nucleon) 


c 


Velocity of the light 


G 


Gravitational constant 



TABLE I: Meaning of the symbols 



We observe that the Gompertz model includes in a very 
interesting and peculiar way nonlinearity, which is a macro- 
scopic mirror of a nonlinear microscopic background ruling 
the growth of natural systems. In fact, the equation looks 
linear in a suitable function y{s) of the original (relative) 
size s, while the nonlinearity is introduced by the fact that 
the function y is itself a nonlinear function of s; in partic- 
ular, y{s) — Ins, or, equivalently, s(y) = expy. On the 
other hand, while the Gompertz equation describes a con- 
tinuous growth of a size in time, we know that many sys- 
tems are placed on some discrete sequence of sizes, ranging 
from a minimum to a maximum scale. Therefore, follow- 
ing the Gompertz model, we try to describe these sequence 
of scales by a sequence {s„} of relative sizes resulting from 
the solution of a linear map fulfilled by a nonlinear function, 
y{sn) = yn, of s„: 
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Vn+l = Snyn, 



(4) 
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where Sn are proportionality coefficients. The solution of the 
map is obviously; 

n 

Vn^yoYi^k- (5) 

Here, s„ will be defined by s„ = z„/zoo, where z„ is the n-th 
size, and z^o is an asymptotic, or reference, size. Moreover, 
for further generalization, we at first allow the dependence 
of the proportionality constant on n (analogously, we could 
promote a to a suitable function of the time in the Gompertz 
equation). The choice of the nonlinear function y(s„) = yn 
is obviously a crucial question: this choice, in principle, will 
depend on the specific system, and it should be made follow- 
ing physical criteria and phenomenological observations. But 
the Gompertz model suggests that the logarithm can play a 
privileged and somewhat universal role. This suggestion is 
strengthened by the connection among Gompertz model, log- 
normality and stability, as we will briefly discuss in the con- 
clusions Js]]. We then select y„ = lns„ (or s„ = expy„). 
Thus, if 7 is defined as in the Gompertz equation, j/o = In 7, 
and the solution (|5]l gives: 

Sn+l^sJ", (6) 

i. e. 

s„^h)^"^=«^''. (7) 

Furthermore, resorting again to the Gompertz model, we make 
the simplifying choice (5„ = (5 = 1 — a, obtaining 

Vn+i = (1 - a)yn, , (8) 
2/„ = (l-a)"ln7,, (9) 

Sn+l=si'-^\ (10) 

and 



In fact, Eq. dSj looks as the form of the discretized Gompertz 
equation. It is, however, worth to be remarked that the phys- 
ical meaning of the discrete map is different from that of the 
continuous equation; in particular, we are not necessarily au- 
thorized to interpret Eq. (O as an equation in (discretized) 
time. 

We now comment on the behavior of the solutions (fTTT i. We 
see that a monotone sequence (increasing if 7 < 1, decreas- 
ing if 7 > 1) is assured only if ci < 1 (while monotonicity is 
always assured by the Gompertz equation). If, instead, a > 1, 
the map ([8]l, (with solution®) changes its sign at each step, 
and, correspondingly, the values of z„ oscillate, but with two 
further possible behaviors. If a is not an integer number, the 
size of the system approaches the asymptotic value z^o by os- 
cillating in alternating way above and below it, with more and 
more damped oscillations. If a is an integer number, the size 



of the system oscillates indefinitely above and below the refer- 
ence value Zoo, assuming in alternating way the pair of values 
and 7^^^"\ and does not converge to any asymptotic 
limit. Therefore, we see that the discrete maps so introduced 
describe a variety of situations. 

In order to investigate the usefulness of the maps ( fTOl i. ( fTTT i. 
we move towards an application, by premising some consid- 
erations. We observe that, if the growth is monotone {a < 1), 
usually one knows (phenomenologically) the initial and the 
final sizes (zo,Zoo) of a system, and that, if 7 < 1 (aggre- 
gation process) the initial value zo is noting but the size of 
the more elementary constituent of the system, and Zqo is the 
maximum size, while if 7 > 1 (fragmentation process), Zoo is 
the size of the more elementary constituent, and zo is instead 
the maximum size. The (phenomenological) knowledge of the 
extremal sizes, together with a good fit for the value of a, is 
sufficient to determine all the intermediate scales. Obviously, 
a more ambitious goal could be to obtain the hierarchical se- 
quence of characteristic scales z„ by fixing phenomenologi- 
cally only the initial size zq, and obtaining the values of a and 
of Zoo by independent theoretical hypotheses. This, as we will 
show, can be done when the physical background underlying 
the system is well established. 

Now, we show that Eq. ( fTTT ) can be applied in the frame- 
work of astrophysics: we assume that the components on dif- 
ferent scales of the observed universe organize themselves ac- 
cording to our scheme, and proceed to verify this hypothesis. 
We choose for z„ the sequence of the length scales ("radii") 
of the astrophysical aggregations, and replace the symbols z„ 
with Rn and Zoo with R, and all the others in a consistent 
manner. Moreover, we fix a = 1/2. Being a < 1, the map 
is monotone. We assume also that Rq is the size of the ele- 
mentary constituent (i. e. the minimum size). Thus, the radii 
Rn are monotonically increasing, and their asymptotic limit 
R is the maximum radius (i. e. the observed radius of the uni- 
verse). We remark that in our computation we consider only 
the order of magnitude of the radii, expressed (in centime- 
ters) as powers of 10. It is natural to choose as elementary 
constituent a nucleon, because the system is ruled by gravi- 
tation, which in turn is determined by the mass distribution; 
and nucleons contain all the significant (observed) mass. Our 
minimum size is then the radius of a nucleon as measured by 
a-particles scattering: Rq = X = 10^^'^ cm. It is also known 
that R, the observed radius of the universe, is i? = 10^® cm 
[4]. Then 7 = lO'^^. The map with 5 = 1/2, becomes 
a square-root map, while if we use s„ = Rn/R we obtain for 
the radii a geometric-mean map: 

i?„+i = (i?„ (12) 

Now, exploiting Eq. ( fTIT i with 7 = lO^''^ and a = 1/2, 
or, equivalently, Eq. ( fT2] i. we obtain the rapidly convergent 
sequence in the first column of TABLE II. 

These length scales represent, in order of magnitude, just 
the main six observed astrophysical length scales ^ (the pres- 
ence of two aggregates both for radius i?i and for radius R2 
is discussed later). Therefore, we have obtained a first im- 
portant result: we can fit the intermediate astrophysical length 
scales by choosing a specific value of a, and by exploiting the 
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RADIUS 


MASS 


AGGREGATE 


10® cm 




planetesimal 


10'' cm 




neutron star 


lO^'' cm 




solar system 


10^** cm 


10«5 


s.m. black hole 


10^^ cm 


lO^^g 


typical galaxy 


lO^^cm 


10«5 


cluster of galaxies 


10^^ cm 


lO^^g 


supercluster of galaxies 


lO^'' cm 


iQ^'g 


observed universe 



TABLE II: List of astrophysical length scales and of the correspond- 
ing mass scales, with the associated astrophysical aggregates 



experimentally determined values of the two extremal length 
scales. 

However, as remarked previously, a scientific theory with 
some predictive power should be able to determine, on the ba- 
sis of few other assumptions, both the (apparently arbitrary) 
choice a = 1/2, and the value of the maximum size R, by 
the mere phenomenological knowledge of the minimum scale. 
We now show that this is actually possible by resorting to 
a suitable assumption naturally suggested by the physics of 
the gravitational systems: We assume that the dimension of 
the aggregate with radius Ri (and mass Mj is given by the 
Schwartzschild radius of a three-dimensional close packing 
of elementary constituents (nucleons: radius X, mass m), i. 
e. by a "collaps" restricted by the constraint that the escape 
velocity is the maximum one: the velocity c of the light. This 
hypothesis leads to the following two conditions: 



M _ m M 



(13) 



where G is the gravitational constant. The second equation 
is the Schwartzschild condition for the radius, while the first 
one requires that the mass density per unit volume of the first 
aggregate coincide with that of a nucleon (three-dimensional 
close packing). By eUminating the mass M, we obtain: 



Ri = (Ai?)5, 



(14) 



with R = {Xc)^/Gm. Inserting the numerical values of 
A, c, g, m provides R = lO^^cm, i. e. just the observed radius 
of the universe. Thus, we see that our simple hypothesis leads 
both to an independent determination of the maximum size R, 
and to the value of a which coincides with our original choice 
(Note in fact that Eq. (fT4l) is Eq. (fT2] | with n = 0). Here we 
comment also that, by defining (in order of magnitude) the age 
T of the observed universe by i? = cT, from the expression 
soon established for R we get: T = X^c/Gm = lO^^s, i. e. 
just the right order of magnitude. 

Finally, we aim to determine, besides the length scales, also 
the masses of the astrophysical aggregates. To this purpose, 
we note that, if some other quantity of the system, say A/„, 
is connected to z„ by the "allometric" relation: Af„ = foz*, 
then the quantity qn = M„/il/oo satisfies the same map as 



s„, and has the same form of solution with the replacement 

7 — >■ 7 = 7'': 



(7 



7 



<5(1-q)" 



(15) 



Moreover, putting together the allometric relations at a generic 
n and at n = 0, we obtain also the proportionality constant b 
as: b ^ Mq/z^. 

In the case of the universe, we consider the sequence of 
the masses A/„, and the sequence of their radii i?„, and find 
their "allometric" relation by introducing another assumption 
(of "minimum fluctuation"), probably less intuitive with re- 
spect to the first one, but reasonable if one considers that the 
very complex character of the gravitational systems, and the 
huge numbers of elementary constituents induce instabilities 
and fluctuations: An aggregate can bind a test particle up to a 
distance where the gravitational force generated by the aggre- 
gate on a test particle reduces to a value comparable with the 
background fluctuating force, whose magnitude is determined 
by the elementary constituent. The claim gives, for each n: 



1. e.: 



bRl 



where 6 = m/X^ = lOOg/cm^. Then, 
(5 = 1/2, we obtain the sequence: 



S = 



Mr. 



7 



2(2" 



-1), 



(16) 
(17) 

2, and, being 
(18) 



where, by denoting M (instead of Moo) the asymptotic value 
providing the total mass of the observed universe, we have 
exploited: m/M = 7^. Then, once the sequence for the 
radii is established, the allometric relation, and the knowl- 
ege of the minimum mass m, automatically provide the se- 
quence of the intermediate mass scales until the maximum 
one. From 7 = 10^"^^ and m = 10"^''^ g, we obtain (in order 
of magnitude) the sequence of the second column in TABLE 
II, which represents just the sequence of the observed astro- 
physical masses. In the last column of TABLE II we describe 
the corresponding structures. Note also that the total mass of 
the observed universe corresponds to that of 7"^ = 10^* nu- 
cleons, in perfect agreement with the central value obtained 
from nucleosynthesis calculations 

We conclude this letter with the following observations and 
remarks: 

Being the macroscopic behavior of complex systems the re- 
sult of the collective behavior of a huge number of elementary 
constituents, it must be a mean effect deriving by a suitable 
microscopic (stochastic) model. In fact, in ref. Isf] it has been 
recently proved that the Gompertz equation is a macroscopic 
consequence of the action of a microscopic, lognormally dis- 
tributed diffusion process performed by the elementary con- 
stituents of the system, because this equation holds for the 
median of the process, which then provides the evolution in 
time of the macroscopic characteristic size. 
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We remark that hypotheses and resuhs on the astrophysi- 
cal structures here discussed, and the hierarchical maps for 
the sequence of astrophysical radii and masses are contained 
in papers previously published by the authors in collaboration 
with other researchers jst]. However, we as well remark that 
in this letter we frame these results in a model based on dis- 
crete maps deduced by the Gompertz model of growth. The 
background stochastic origin of the Gompertz model [3l, to- 
gether with the derivation from physical hypothesis of the val- 
ues of a and of the asymptotic length R, lead us to interpret 
our scheme as a first important step towards the discovery 
of the underlying process performed by the elementary con- 
stituents, and responsible for the observed sequence of astro- 
physical scales. Note, in fact, that it clearly emerges from our 
model the role of a geometrical progression, whose relevance 
in the framework of natural phenomena and in a variety of ex- 
perimental environments was highlighted already at the end of 
nineteenth century by F. Galton and by D. McAlister 
which showed that the geometrical mean (median) describes 
the behavior of a large set of natural phenomena better than 
the arithmetic one. 

The discrete maps here introduced can provide a general 
model that can be potentially exploited to search for the en- 
lightenment of a large number of growth phenomena in many 
fields of research. Namely searching for sequences of sig- 
nificant scales associated to different systems, which possi- 
bly develop these intermediate scales step by step in time, or 
which organize themselves on all these scales simultaneously. 
It is also of great interest to investigate the existence of natural 



quantities whose evolution is characterized by oscillations, as 
described by the maps with a > 1. 

We remark that the apparent discrepancy between the first 
relation in Eq. (fTsT i and the relation (fT6] l is solved by the fact 
that to the same radius can be associated the mass Af„ and 
also the mass Mn+i, where the mass i\/„ satisfies the relation 
(fT6] l (two-dimensional close packing), while the mass Af„+i 
satisfies the relation (fOT i (three-dimensional close packing, 
"critical" aggregate) ISfl. For example, as is well known, 
a neutron star (three dimensional close packing with mass 
M2 = lO'^'^g) has, in order of magnitude, the same radius 
[Ri = lO^m) of a planetesimal (two dimensional close pack- 
ing with mass Mi = lO^^.g). Similar considerations can be 
extended by comparing a typical galaxy with a supermassive 
black-hole at the center of the galaxy bulge, whose estimated 
mass (until billions of solar masses, i. e. lO^^g) is in fact com- 
parable with that of the whole galaxy This explains the 
two structures associated to the first two radii in Table I. Sum- 
ming up, our model allows two possible sequences, given by: 
{Rn,Mn} (two-dimensional close-packed aggregates) and 
{Rn, Mn+i} (three-dimensional close-packed, critical aggre- 
gates), and this scheme is confirmed at least on scale and 
on scale i?2. 

A question which is worth to be deepened is the meaning of 
the proportionality parameters in the linearized maps. But this 
question is connected to the more general, and already raised, 
question of the possible stochastic background underlying the 
maps, and will be addressed in forthcoming papers. 
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